KRALL-JACOBI COMMUTATIVE ALGEBRAS OF PARTIAL 
DIFFERENTIAL OPERATORS 



PLAMEN ILIEV 



Abstract. We construct a large family of commutative algebras of partial 
differential operators invariant under rotations. These algebras are isomor- 
phic extensions of the algebras of ordinary differential operators introduced 
by Griinbaum and Yakimov corresponding to Darboux transformations at one 
end of the spectrum of the recurrence operator for the Jacobi polynomials. The 
construction is based on a new proof of their results which leads to a more de- 
tailed description of the one-dimensional theory. In particular, our approach 
establishes a conjecture by Haine concerning the explicit characterization of 
the Krall- Jacobi algebras of ordinary differential operators. 
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1. Introduction 

It is practically impossible to list the numerous applications of Jacobi polyno- 
mials which were introduced more than 150 years ago as solutions of the hypergeo- 
metric equation [10]. In 1938, Krall [13] studied the general problem of classifying 
orthogonal polynomials which are eigenfunctions of a higher-order differential oper- 
ator and solved it completely for operators of order 4, thus extending the classical 



Date: March 21, 2011. 

The author was supported in part by NSF Grant #0901092 and Max-Planck-Institut fur Math- 
cmatik, Bonn. 

1 



2 



P. ILIEV 



orthogonal polynomials [14] , In the last century, many different solutions of Krall's 
problem were constructed, see for instance [IT] [121 HH US] and the references therein. 

More recently, new techniques have emerged in the literature inspired by the 
bispectral problem 2 ]. The most general result concerning the Krall's problem 
appeared not long ago in the beautiful work of Griinbaum and Yakimov [5]. They 
constructed a large family of solutions to a differential-difference bispectral prob- 
lem, containing as special or limiting cases all previously known solutions of Krall's 
problem. Their approach is based on a very subtle application of the Darboux 
transformation [JJ, one of the basic tools in the theory of solitons PJj]. The con- 
struction naturally leads to a commutative algebra of ordinary differential operators 
diagonalized by the generalized Jacobi polynomials, which we call the Krall-Jacobi 
algebra. 

In the present paper we discuss an extension of the above theory within the 
context of partial differential operators invariant under rotations. First, we give a 
new proof of the results in [5] , which shows that if we iterate the Darboux trans- 
formation at one end of the spectrum of the Jacobi recurrence operator, then the 
Krall-Jacobi commutative algebra of differential operators is contained in an as- 
sociative algebra with two generators which have natural multivariate extensions. 
This fact does not seem to follow easily even from the explicit formulas in [TT] for 
the operator of minimal order in the Krall-Jacobi algebra in the case of a single Dar- 
boux transformation. Our proof is based on the approach used in [8] [9] to establish 
the bispectrality for rank-one commutative algebras of difference or q-difference op- 
erators. The main difficulty here is to evaluate certain discrete integrals (or sums) 
involving the Jacobi polynomials, which were trivial integrals involving exponents 
in the rank-one case. As another corollary of our construction we establish Con- 
jecture 3.2 on page 161 in [6] for the Krall-Jacobi algebra of ordinary differential 
operators, which gives an explicit characterization of the isomorphic (dual) algebra 
of eigenfunctions. Moreover, our techniques yield an explicit eigenbasis of multi- 
variate polynomials for the Krall-Jacobi algebras of partial differential operators in 
terms of the quantities used to describe the sequence of one-dimensional Darboux 
transformations and the spherical harmonics. 

The paper is organized as follows. In Section [2] we introduce the Jacobi polyno- 
mials and the corresponding recurrence and differential operators. In Section [3] we 
review briefly the sequence of Darboux transformations from the Jacobi recurrence 
operator. In Section [5] we present the new proof of the results in [5] together with 
the additional properties of the Krall-Jacobi algebra mentioned above. We give 
a detailed proof for the case needed for the multivariate extension (i.e. Darboux 
transformations at one end of the spectrum), but we indicate the necessary mod- 
ifications for Darboux transformations at both ends in Remark 14.61 In Section [5] 
we define the Krall-Jacobi algebras of partial differential operators and we write an 
explicit basis in the space of polynomials in several variables. In the last section we 
illustrate the constructions in the paper with the simplest possible example which 
leads to a multivariate analog of Krall polynomials [14] . 
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2. JACOBI POLYNOMIALS AND OPERATORS 



Throughout the paper we denote by p" (z) the Jacobi polynomials normalized 
as follows 



Pn' P (z) = (-1) 



(q + /3 + l)n F 



(-n,n + a + + l,0 + l;z), 



(2.1) 



where F is the Gauss' 2-F1 hypergeometric function and (a) n is the the shifted 
factorial (a)o = 1, (d) n = a(a + l) • • • (a+n— 1) for n > 0. Besides the orthogonality 
on [0, 1] with respect to the measure (1 — z) a z^dz the Jacobi polynomials are 
eigenfunctions of the second-order differential operator 

B a ^(z,d z )=z(z-l)d 2 z + (z(a + f3 + 2)-((3 + l))d z (2.2) 

with eigenvalue 

X" +l3 = n(n + a + + 1), (2.3) 

i.e. we have 

B ai0 (z,d z )p^(z)=X^p^(z). (2.4) 
In order to simplify the notation, we shall write A n instead of A" +/3 unless we need 
to use different values for the parameters and the explicit dependence on a + is 
important. 

Since any family of orthogonal polynomials satisfies a three-term recurrence re- 
lation, the polynomials p^{z) are also eigenfunctions for a second-order difference 
operator acting on the discrete variable n. More precisely, if we denote by E n the 
customary shift operator acting on functions of a discrete variable n by 

E n fn fn+Xi 

then 

L a ^(n,E n )p^(z) = zp^(z), (2.5) 
where L a ^(n, E n ) is the second-order difference operator 

L a ,p(n t E n) = A n E n + B n ld + C n E~\ (2.6) 



with coefficients 



C, 



(n + l)(n + /3 + l) 



(2n + a + (3 + l)(2n + a + (3 + 2) 
(n + a)(n + a + 0) 



(2.7a) 



(2.7b) 
(2.7c) 



(2n + a + 0) (2n + a + /3 + 1) 
B n = 1 — A n — C n . 

Here and later we assume that p"'^(z) = for n < 0. Equivalently, if we think of 
p^{z) as the semi-infinite vector 

\p^( Z ),p^(z),p^(z),...r, 

then L ai p(n, E n ) can be represented by the tridiagonal semi-infinite (Jacobi) matrix 

'B A 
Ci B x A± 

L a , p (n,E n )= c 2 B 2 A 2 ■ (2-8) 



In view of equations (|2.4|) and (|2.5j) we can say that the polynomials p"'°(z) solve 
a differential- difference bispectral problem. 
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Using the explicit formula (|2.1[) one can easily check that the Jacobi polynomials 
satisfy the following differential-difference equation 

[2zd z + (a + P)m^z)-p^M) = (2n + a + /3)(p?/ (z) + p^z)). (2.9) 

The above formula will be a key ingredient in the computation of "discrete integrals" 
involving the Jacobi polynomials, which are analogous to standard (continuous) 
integrals of products of polynomials and exponents. This is the first place where 
the particular normalization of p"' (z) is very important. 

3. Discrete Darboux transformations 

In this section, following [5], we describe the result of k successive Darboux 
transformations starting from L a ^(n, E n ) at z = 1 for a £ N and k < a. 

Recall that the lattice version of the Darboux transformation 16J of a difference 
operator (bi-infinite matrix) Co at z amounts to performing an upper-lower factor- 
ization of Cq — zld and to producing a new matrix by exchanging the factors. If 
we iterate this process k times starting from Cq at z = 1, we obtain a difference 
operator C as follows: 

Co = Id + VoQo r\ d = Id + Q a Po = Id + ViQi rv - • • 

£ fc _i = Id + Q k - 2 Pk-2 = Id + Pk-iQk-i (3-1) 

r\C = C k =Id+ Qk-iVk-i- 

From ([3~Tj) it follows that 

CQ = QC (3.2) 

and 

(Co - W) fc = VQ, (3.3) 

where V = VoV\ ■ ■ - Vk-i and Q = Qk-iQk-2 ■ • • So- Tbe above formulas imply 
that 

ker(Q) C ker((£ - W) fe ) and £ (ker(Q)) C ker(Q). (3.4) 
Conversely, one can show that if (|3.4[) holds then there exists a difference operator 
(bi-infinite matrix) C such that (|3.2p holds and this operator is obtained by a 
sequence of Darboux transformations as in (|3.1[) ■ 

Note that, up to a factor, the lower-triangular bi-infinite matrix Q is uniquely 
determined by its kernel, i.e. if {V*n^}*=o is a basis for ker(Q) then 

Qfn = .g n Wr„(^ 0) , ^\ . . . , 4 fe_1) » fn), 

for an appropriate function (diagonal bi-infinite matrix) g n . We use Wr„ to denote 
the discrete Wronskian (or Casorati determinant): 

Wv n (h£\h£\ h^) = det(/i» whKijKk. 

Combining the above remarks, we see that the sequence of Darboux transformations 
(|3.ip for Co is characterized by choosing a basis for ker(Q) satisfying 

(Co ~ Id)^ 0) = and (C - Id)^ = for j = 1, . . . , fc - 1. (3.5) 

We apply the above construction for the difference operator Co = C a ,fs(n + 
e, E n ) = A n+e E n + B n+5 ld + Cn+eE" 1 . Note that for a generic e, the coefficients 
A n+£ , B n+£ and C n + e are defined for all n £ Z and therefore C a ,p(n + e, £"„) is a 
well defined second-order difference operator (or a bi-infinite tridiagonal matrix) . 
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For j S {0, 1, . . . , fc — 1} we set 

— n — n/ — R \ i 

(3.6a) 



ld _ (-lp(n + l) J (-n-g-/3) J 



,2,j _ + l)g(n + (3 + l) a (-n),-(rc + a + /? + 1), 

" " ]\a\{l + a) 3 {l + P) a ■ (3 ' 6b) 

One can show that the functions 4>n J +e are linearly independent (as functions of n) 
and 

(£ a ,p(n + e, E n ) - Id)#£ £ - 0, for i = 1,2 

(n + £,£„) -Id)<^_ £ = ^e 1 ) for j = l,...,fc- 1, £=1,2. 

Thus, we can write the functions ^« as a linear combination of <filf +e as follows 

i 

} = ^(aj-l^n+e + bj-l'Pn+e)- 
1=0 

If bo = then ipn^ = clq ^ 0, i.e. we can take f/'i ' = 1 and one can check that 
C\ in (|3.ip coincides (up to a conjugation by a diagonal matrix) with L a -\^(n + 
e,E n ). Therefore the operator C — Ck can be obtained by a sequence of fc — 1 
Darboux transformations starting from C a -i,@{n + e,E n ). Thus we can assume 
that bo 0, hence we can take bo = 1. Since Q depends only on the space 
span{i/)i°\ ipn \ ■ ■ • , ipn }, and not on the choice of the specific basis, we can can 
take bj = for j > 0. Thus we shall consider a basis for ker(Q) of the form 

3 



yv ■'--> 

1=0 

depending on a e N, /3 and fc free parameters ao> a%, . . . , a^-i- 

Now we consider the limit e — > 0. Since lim e ^o = 0, it follows that the 

intertwining relation (|3.2j) holds for the semi- infinite parts of the bi-infinite matrices 
Co, C and Q. In other words, we have 

L{n, E n )Q = QL a< p(n, E n ), (3.7) 

where L a ^(n, E n ) is the semi-infinite matrix in (|2.8p . L(n,E n ) is a similar semi- 
infinite Jacobi matrix, and Q is a lower-triangular semi-infinite matrix, acting on 
functions (vectors) /„ = [/ , /i, . . .]* by 

Q/„ = g„Wr„(Vi 0) , ■ • .,^* -1) ,/n), (3-8) 

where 

^^^a^^ + ^i, (3.9) 
;=o 

with the convention that /j = when j < 0. 

We shall normalize the matrix Q by taking g n — 1 in (|3.8I) . Finally, we denote 
by q n {z) the polynomials defined by 

q„(z) = Qp^(z) = Wr n (^\^\...,^ k - 1 \p^(z)) ) (3.10) 

which depend on the free parameters a <E N, (3 and a = (ao, ai, . . . , afe-i). 
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Remark 3.1. To simplify the notation we omit the explicit dependence of the 
parameters a, f3 and a in the functions ipn^ and q n (z). When these parameters are 
needed, we shall write ipn " > ' a '^' a [z) and q" ,l3 ' a (z). 

Remark 3.2. From (|23]) and (|3~7| it follows that 

L(n, E n )q n (z) = zq n {z). (3-11) 

It is easy to see that the off-diagonal entries of the tridiagonal matrix L(n, E n ) are 
nonzero and therefore, by Favard's theorem, there exists a unique (up to a mul- 
tiplicative constant) moment functional A4 for which {<7n(z)}^Lo i s an orthogonal 
sequence, i.e. 

M(q n q m ) = 0, for n ^ m and M(qfc) ± 0. 

More precisely, one can show that there exist constants Uo,u\, . . . , itfc_ i such that 
the moment functional A4 is given by the weight distribution 

fc-i 

W( Z ) = (1- z) a - k Z P +J2 U 3 SU) (Z -1), 

where <5 is the Dirac delta function. The parameters Uj correspond to a different 
parametrization of the Darboux transformation (|3.ip . The proof of this statement 
is well known and we omit the details. We refer the reader to jU Theorem 2, p. 287] 
where a similar result was proved for extensions of Laguerre polynomials. 

4. The commutative algebras A a ^' a and T> a - j3 ' a 

In this section we prove that the polynomials q n {z) are eigenfunctions for all 
operators in a commutative algebra of differential operators. 

4.1. Statement of the main one-dimensional theorem. In order to formulate 
the precise statement let us denote 

£>! =zd z &ndV 2 = zd 2 z + {(3 + l)d Zl (4.1) 

and let 55 ^ denote the associative algebra generated by T>\ and T> 2 , i.e. 

J)p=R(V u V 2 ). (4.2) 

It is easy to see that 

[V 2 ,V 1 ]=T>2. (4.3) 

Next we set 

r n =Wr n (tf)tf,... i e i) ) (4-4) 

where the functions ipn ^ are given in p.9[) . Here we use the same convention as in 
Remark l3.1l and we omit the explicit dependence of the parameters a, j3 and a. Note 
that T n defined above and A„ defined in (|2.3I) are polynomials of n. In fact, one can 
show that up to a simple factor, r n belongs to R[A"^ fe _ 1 ^ 2 ] = R[A" +,9 ~ fc+1 ], where 
A* = n(n + s + 1). One simple way to see this is to use the involution introduced 
in [S] which characterizes the subring R[A^] in M[n]. For s e R we define /( s ) on 
R[n] by 

I {s) (n) = -(n + s + 1). 

Then clearly I {s) (X s n ) = X s n hence every polynomial of A* is invariant under the 
action of I*- 5 - 1 . Conversely, if r G R[n] is invariant under I^ s \ then r g K[A^]. 
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Note also that the functions are invariant under l( a+ P\ and therefore ipn^ 
in (|3.9|) are also invariant under /( Q +^). From this it follows that ji a +P~ k + 1 ) w [\\ 
reverse the order of the rows in the determinant in equation (|4.4j) . leading to 

= ( _ 1 jfc(fc-l)/2 7H- 

The last formula shows that if k = 0, 1 mod 4, then r„ £ R[A„_(/ c _i)/2]- Otherwise, 
t„ is divisible by (2n + a + /3 — fc + 2) = \ n -k/2+i ~\i-k/2 an d the quotient belongs 
to K[A n _(fc_!)/ 2 ]- Summarizing these observations we see that 

T n = el fc) T = (A„_(fc-l)/2), 

where f is a polynomial and 

4fc) = jl if A; = 0,1 mod 4 

[A n _ fc / 2+ i - A n _ fc /2 iffc = 2,3 mod 4. 

For n, m £ Z and for a function / s defined on Z it will be convenient to use the 
following notation 

f s d^ d (s) = { if n = m 

E^Ln+i fs if " < m. 



Thus 



/>n />n— 1 

/ f s dn<i{s) = fn+ / fsdfid(s) for all m, n £ Z. 



(4.6) 



We denote by .4 a '^ a the algebra of all polynomials / such that /(A n _fc/ 2 ) 
/(A„- fe /2-i) is divisible by t„_i in R[n]: 



^.fta = g R[t] . j(W /(*n-fc/2-l) g R[ 



T n -1 



(4.7) 



Remark 4.1. It is not hard to see that A a '^' a contains a polynomial of every 
degree greater than deg(f). Indeed, notice that for every r £ R[f] there exists 
f £ l[t] such that 

r(A n ) - r(A»_i) 

t r = J"(A„-i/ 2 ). 

Conversely, for every polynomial f there exists a polynomial r such that the above 
equation holds. Moreover, up to an additive constant, r is uniquely determined by 

pn pri—1 

r(X n ) = J (A 8 - A a _i)f(A a _i/2)d/id(s) = y (A s+i - A s )f(A s+ i/ 2 )rf^ d (s). 

If f ^ then deg(r) = deg(f) + 1. Thus, for every g £ R[t] and c £ R we can define 
/ £ A a ^ a by 

/(A„_fc/ 2 ) = / ei fc+2) g(A s _ (fe+1)/2 )r s _id/i d (s) + c. 
Jo 

Conversely, to every / £ A a '^' a there correspond unique g £ R[t] and c £ R, so 
that the above equation holds. 

The main result in this section is the following theorem. 
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Theorem 4.2. For every f 6 A a ^' a there exists Bf = BjiVi^V-i) G £>/3 such that 

B f q n {z) = /(A n _ fe/2 )g„0). (4.8) 

Thus, V a -P'- a = {B f : / € ^"• /3;a } zs a commutative subalgebra ofQp, isomorphic 
toA a ' f3 ' a . 

Equations (|3.11[) and (|4.8|) establish the bispectral properties of the generalized 
Jacobi polynomials q n {z). The fact that for every / £ _4 Q ' ,9;a there exists a differ- 
ential operator Bf satisfying (|4.8p was conjectured in [BJ. The fact that Bf G S,g 
is crucial for the multivariate extensions. 



4.2. Auxiliary facts. For the proof of the above theorem we shall need several 
lemmas. First we formulate a discrete analog of a lemma due to Reach [T7] . 

Lemma 4.3. Let /n ,/n ,•-•)/« oe functions of a discrete variable n. Fix 
ni, ri2, • • • , Uk+i S Z and Zei 



k+i 



^=E(-i) w+j /« / /rwr s (/^,...,/w,...,/r +ij )^ d (s), 



t(i) 



p(j) 



f (fc+i)> 



with the usual convention that the terms with hats are omitted. Then 
Wr„ (/«, . . . , fW,F n ) = / /^Wr^/P, . . . , /W)d W (a) 

xWr^/W,...,/^). 



(4.10) 



The above lemma was used in [8] to give an alternative proof of a theorem in 
[7] that rank-one commutative rings of difference operators with unicursal spectral 
curves are bispectral. Similar argument was used also in |9 for rank-one commuta- 
tive rings of g-difference operators. Since the application in our case is very subtle 
and we need different elements of the proof, we briefly sketch it below. 

Proof of Lemma \4-3\ Note that 



0, for every I = 0, 1, . . . , k — 1. 



Expanding the above determinant along the last row we obtain 

fe+i 

^ ( _ 1)fe+ i +J - / (i)^ Wrn(/ (i )j _ _ ; f (k + i) } = o for / = 0, . . . , A - 1. (4.11) 

Using (|4.6I) and (14. 1 1[) we see that for I = 0, . . . , k — 1 we have 

fc+l r. n 

F n . l = ^(-l) k ^f^ l / f^Wr s (fP,...J^\...j( s k ^)df, d ( S ), (4.12) 

.7=1 ^ 



Jn 


f(2) 
,/n 


,(fc+l) 


Jn-1 


f (2) 
./n-1 


,(fc+l) 


,(1) 

Jn-fc+l 


,(2) 

Jn-fc+l 


f (fe+l) 
Jn-fc+l 


f U) 
Jn-l 


f (2) 


„(fc+l) 
■Jn-Z 
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and 

fc+1 ~n 



(4.13) 



/^Wr^/W,...,/^ 1 )). 



If we plug gUg) and ([4TT3"|) in Wi n ($\ . . . , f n k \ F n ), then most of the terms cancel 
by column elimination and we obtain (|4.10l) . □ 



Remark 4.4. We list below important corollaries from the proof of Lemma [ 

(i) Note that the right-hand side of (|4.10p does not depend on the integers riy , n% . . . , n^. 
Moreover, if we change n^+i then only the value of 

rn-l 

/i )Wr s (/W,...,/i fc ))^ d ( S ) 

will change by an additive constant, which is independent of n and f n k+1 ^- Thus, 
instead of (|4.9p we can write 

fc+1 r-n 

F n = J2(- 1 ) k+1+i fi j) / /^Wr^/W,...,/^,...,/^ 1 ))^^), 

.7=1 J 

leaving the lower bounds of the integrals (sums) blank and we can fix them at 
the end appropriately. This would allow us to easily change the variable, without 
keeping track of the lower end. 

(ii) From (|4.11l) it follows that for every I = —1,0,1 ... ,k — 1 we can write F n also 
as 

fc+1 pn+l 

F n = / /i 0) Wr s (/«, ...,%),..., fW)dto(s), 

3=1 J 

and changing the variable in the discrete integral we obtain 

fc+1 f*n 

F n = £(-l) fc+1+ V« / /fflWr s (/« , . . . , f% . . . , / s ( + ! 1) )^d( S ). 

3=1 

(iii) Suppose now that k is even. Using (ii) with / = k/2 — 1 we have 
fc+i 



F - Vf-l) H1+3 fW / f (0) Wr (f (1) f U) f (fe+1) )duA(s) 

rn — L > Jn / Js+k/2-1 VVls Us+fc/2-l' ■ • ■ ' Js+fc/2-1' ■ • ■ ' J s+k/2-1 ) a ^ d \ b ) ' 

(4.14) 



3=1 



Let us consider the sum consisting of the first k integrals: 

fc />n 

p(fc) = Vf-11 fe+1+J ' f W / f (0) Wr (f {1) f U) f (fe+1) )duj(s) 

3=1 

Expanding each Wronskian determinant along the last column we can write Fn 
as a sum of k terms F n k,m \ each one involving as integrand one of the functions 
/i+ra j where m = —k/2, —k/2 + 1, . . . , k/2 — 1. We can use (|4.11[) once again, 

this time for the functions f n °\ fn \ • • • , fn (i-e. omitting f n k+1 ^), to change s as 
follows: 
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• If to > we can replace n with n — m in the upper limit of the integral, 
or equivalently, if we keep the upper limit of the integral to be n, we can 
replace s by s — m in the integrand. Thus F^ k,m ^ will have /i fc+1 ' as 
integrand (and the integration goes up to n) 



• If m < —1 we can replace s by s — to — 1, thus F)t'" l> will have /j*^ 
integrand (and the integration goes up to n). 

This means that we can rewrite Fn as sums of integrals, and the integrands can be 
combined in pairs (corresponding to m and (— to— 1) for to = 0, 1, ... , fc/2) involving 



and /„ 



(fc+i) 



Explicitly, we can write Fn 1 , as a sum of terms which, up to a 



sign, have the form 



J n 



(i) 



.(0) 

s-m+k/2-1 



f(2) 

s+fe/2-m-l 
f(2) 

s+fc/2-m-2 



f(2) 



f (2) 

■* s — k/2—m 
(2) 



r(*0 

s+fe/2-m-l 
(fc) 

s+fc/2-m— 2 



i(k) 



f. 



(0) 

s+m+k/2-1 



s+fe/2+m-l 
,(2) 

s+fc/2+m— 2 



m 

Js-1 



s—k/2—m 

,(ao 



P (fc+i) 



s+fc/2+rn-l 
•(fc) 

s+fc/2+m-2 



f(fc) 
/s-1 



(2) 

s — fc/2+m 



(fc) 

S-fe/2 + ?Tl 



,(fc+l) 

/s-1 



d/i d (s), 



(4.15) 



for to = 0, 1, ... , fc/2. For simplicity, we wrote explicitly only the terms with f$p in 
front of the integral, but we have also similar expressions obtained by exchanging 
the roles of f n ^ and fn^ for every j = 2, . . . , fc. 

(iv) Finally if fc is odd, we shall use (ii) with I = (k — l)/2 and I = (k — 3)/2 and 
write F n as the average of these two sums. Then we can apply the same procedure 
that we used in (iii) to write F^' as sums of integrals, whose integrands can be 
combined in pairs involving fi k+1 ^ and / s _^ . 



The second important ingredient needed for the proof of Theorem 
following lemma. 

Lemma 4.5. For r(n) £ M[n] there exists B £ T)^ such that 

n 

[r(s)p^(z) ~ r(-s ~a- (3)p a s L\(z)W A (s) = Bp^(z). 

-l 

Proof. Let us rewrite r(n) as a polynomial of 2n + a + /3, i.e. we set 



the 



(4.16) 



r{n) = r(2n + a + /3). 
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Then r(— n — a — (3) — r(—(2n + a + /?)). Thus, it is enough to show that for every 
j G No there is Bj G Dp such that 

(2n + a + PY (p%P(z) + (-ly+'p^iz)) = Bj (p^(z) - p^z)) . (4.17) 

Indeed, if (|4~T71) holds then 

fjls + a + (3Y (pf'\z) + (-ly+^iz)) dfx d (s) 

/n 
{pf( 3 ( Z )-p a s f l{z)Wd{s) 

giving the proof when f{t) — P , hence for arbitrary polynomials by linearity. When 
j is odd, (|4.17j) follows immediately from (|2.9[) and the definition of Dp, see equa- 
tions (|4.ip ~ (|4.2p . The case j — is obvious and therefore it remains to prove the 
statement when j > is even. Note that (2n + a + /3) = A„-A n _i. Thus (|4T7l) for 
j even will follow if we can show that there exist operators B',B" G Dp such that 

(A„ + An_i) (p^(z) -p:A{z)) = & (p^{z) -K-to) (4.18a) 

and 

A n A„-i (p^(z) -p^z)) = B" (p«;'\z) -p^z)) . (4.18b) 

Using (|4.1I) we see that the Jacobi operator B a> p(z,d z ) defined in (12 .2|) belongs to 
Dp since 

d z ) = V\ + (a + P + 1)2?! - XV (4.19) 
From equations (12 .4[) and (|2.9[) we find 

A„-iK^( z ) - KPn'A(z) =(K-i - K) ( P a J{z) +P^ 1 (Z)) 

+ \ n p^{z)-\ n - l p a J 1 (z) 
=B'"{p^(z)-p^ 1 (z)), 

where B'" = B a ^(z, d z ) — 2£>i — (a + (3) G Dp. Using the above equation together 
with (|2.4p we obtain 

(A„ + A„_0 ( P f/(z) -p£?i(*)) =\ nP ^{z) Xn-xp^iz) 

+ \ n _ 1 p^(z)-\ n p a J 1 (z) 
^B'(p^{z)^p a J 1 {z)) 1 
where B' = B aj/3 (z, d z ) + B'" G Dp proving (|4.18al) . Similarly, 

A„A„_! [p^{z)-p a n ^)) =B a A*,9*) - A„K-V*)) 

^"(K'^W-P^iW), 

where B" — B a ^p(z 1 d z )B"' G 23/3, establishing (|4.18b|) and completing the proof. 

□ 
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4.3. Proof of the main one-dimensional theorem. We are now ready to give 
the proof of Theorem 14.21 Let / 6 A a '^' a . From Remark 14.11 we know that, up to 
an additive constant, we have 

/n— 1 
e i+l 2) .9(^s-(fc-l)/2)T ;5 C?Md(s). 

We apply Lemma 14.31 and Remark 14.41 (i) with 

Jn - e n+l 9{\i-(k-l)/2) 

fi j) =rtr 1] fori = 1,2,...,* 

Using equations (|3 . 10[) and (I4.4[) we see that the right-hand side of (|4.10[) is equal 
to (/(A„_fc/ 2 ) + c)q n (z), where c is a constant (independent of n and z). The goal 
now is to show that if we choose appropriately the integers rij in Lemma 14.31 then 
there exists a differential operator Bf G Dp such that 

F n =B f p a J{z), (4.20) 

Suppose first that k is even and let us write F n as explained in Remark 14.41 (iii) . 
From Remark 14.11 it follows that the integral in the last term in the sum (I4.14[) 
representing F n is an element of M[A n ]. Therefore, using (|2.4I) and (14. 19)) we see 
that the last term in this sum is of the form Bp^{z) for some operator B G Dp. 

(k) 

Thus we can consider the sum Fn of the first k terms. It is enough to show 
that each term of the form (|4.15[) can be represented as Bp"'@(z) for some operator 
B G Dp. Recall that fn G R[A„] for j = 1,2, and therefore it suffices to 

show that the integral is of the form Bp"'"(z) for B G Dp (since then we can 
commute B and /„ and use (|2.4p ). Now we apply Lemma [4.51 To simplify the 
argument, we shall use the involution /( Q +' 3 ~ 1 ) which acts on polynomials in R[n] 
by I^+P-^in) = -(n + a + f3). Thus for r(n) G M[n] we have I i - a+l3 -^{r{n)) = 
r(— (n + a + /?)). Note that for every / G R we have 

/(^- 1 )(A n+i ) = A n _ i _ 1 . 

Therefore, if we denote by det ' n and det n the determinants in f|4.15[) . then 
will reverse the order of the rows, hence 

l( a+ P-V(det' n ) = (-l) (k -^ k ^ 2 det' n . 

Since 

f(0) _ (k+2) n x 

Jn±m+k/2-l ~ t n±m+k/2yy^n±m-l/2)i 

and 

^ (a+/3_1) (g(A n+m -i/2)) = g{Ki- m -i/2), 

it remains to check that 

j[a+P-\)( ,(fe+2) n _ /_ n (fe-l)(fe-2)/2 (fc+2) 

1 \ c n+m+k/2' V L ! c n-m+fc/2' 

which follows at once from (|4.5[) by considering the two possible cases k = mod 4 
and k = 2 mod 4. 

The case when k is odd can be handled in a similar manner, using Remark 14.41 
(iv). □ 
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Remark 4.6. If /3 € N then instead of the Darboux transformations (|3.1[) at z = 1 
we can consider a sequence of Darboux transformations at z — 0. Theorem 14 . 21 can 
directly be applied in this case by replacing z with 1 — z and exchanging the roles of 
a and /3. In particular, the corresponding commutative algebra T> a ^' ,a constructed 
in Theorem 14.21 will be a subalgebra of D a = R(£>i, 2?2), where 

T) 1 = (z- l)d z and 2% = (1 - - (a + 1)<9 Z . 

Finally, if both a and /3 are positive integers then we can iterate the Darboux 
transformation at z = 1 and z = at most a and /3 times, respectively. In this 
case, we need also the functions 

^ = jlaKl-PUn + 1),, ' (4 ' 21a) 

Vn = WHWTWi ' (4 ' 21b) 

which are linearly independent and satisfy 

C at p(n + e,E n )<Pn +e = 0, for i = 1,2 

C a ,p[n + e, E n )ip l n 3 +£ = Vn+e 1 , for i = 1, 2 and 1 < j < 0. 

Note that if we define 

t _, n« ( w + 1 )g 
(n + 1) Q 

then S,n i Pn J become polynomials of A„. If we apply k Darboux steps at 1 and I 
Darboux steps at 0, then instead of (|4.4[) we define r„ by 

r n = (Un + a-k-l + 2) k+l _ 1 ) 1 Wr„(^ 0) , 1>g\ . . , i>t% 

where ipii are defined similarly to ip„ f , using the functions <y9^ J . Then Thcorcm l4.2l 
holds with Dp replaced by D = R(zd z ,d z ), since we need to use now operators from 
both algebras Dp and D a . The proof follows along the same lines, using Lemma 
14.51 together with the analogous statement concerning the transformations at 0, i.e. 
for r(n) £ M[n] there exists B € D a such that 

1 f n 

* ' r r(s)Z s p^(z) - r(-s -a- ^-xp^z^d^s) = Bp a J(z). 



6, 



1 



Note that the construction of the operator Bf(T>\,T>-i) in the proof of Theorem 
|4~2] depends only on equations ([2^4]) . (|2~9l) and (|4.19J> . For s e No let us denote 

and let us define B a ^p_ s {z, d z ) similarly to B ay p(z, d z ) using (|4. 19|) with T>2 replaced 
by T> 2 ,s- 

B a ,p, s (z, d z ) = V\ + {a + p + 1)2?! - V 2tS . (4.23) 
Then it is easy to check that 

z- s ' 2 B a ^ s {z,d z ) ■ z s ' 2 = B a , p+S (z,d z ) + A«+f /2 - \°+P +s . 

From this relation and (|2.4p it follows immediately that 

B a ^ 8 {z,d z ) \^+'{z)M = \^ s/2 p^ +s (z)z^. (4.24) 
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It is also easy to see that 

(2zd z + a + f3)(p^+ s (z)z^+p^+ s (z)z^) 

= (2(n + s/2) + a + f3)(p^ +s (z)z s / 2 -p a J+ s (z)z^). 
Note also that 

[Pa, a ,2?i] =V 2 , a , 

which shows that for every / G A 01 ' 13 ' 0, we have a well-defined operator Bf(T)i,T>2, s )- 
Comparing equations f[2T4|) . ((2T9J) and (j47T9|) with (|4~24|) . (|425]l and (f4T23|) we obtain 
the following corollary of Theorem 14.21 



Proposition 4.7. // we define for s G No 

= Wr„(^°i/f ■ • ■ , C" s /2^ ; >n (*)), (4-26) 
i/ien /or every f G J[ a '^' a we have 



Bf{p u V 2i „) i$g»{z)z'/*\ =f(Kl(s-k)/2)%f a (z> S/2 , ( 4 - 27 ) 



where Bf are the operators constructed in Theorem^ 

We have displayed all parameters in (|4.26[) to underline the fact that in the 
functions tpn the parameters stay the same, only the variable n is shifted by s/2, 
while in p"'@(z) we change only the parameter 0, 

5. Krall-Jacobi algebras in higher dimension 

In this section we show that the commutative algebra T> a '^' a of ordinary differen- 
tial operators constructed in Theorem 14.21 is isomorphic to a commutative algebra 
of partial differential operators invariant under rotations which can be diagonalized 
in the space of polynomials in d variables. 

5.1. Notations. Let x = (x\,X2, ■ ■ ■ ,Xd) G K d , and let V d = M.[x] = R[xi,X2, ■ ■ ■ ,Xd] 
be the corresponding ring of polynomials in the variables x\, X2, ■ ■ ■ , Xd- We denote 
by B d and S' d_1 the unit ball and the unit sphere in M. d : 

B d = {x G R d : \\x\\ < 1}, S^ 1 = {16^: ||x|| = 1}. 

In polar coordinates we shall write x = px' where p = ||x|| and x' G S^ -1 . 

We denote by A x = J2j=i ®x tne Laplace operator and by Hi the space of 
homogeneous harmonic polynomials of degree I, i.e. the homogeneous polynomials 
Y(x) of degree I, satisfying the equation A x Y(x) — 0. It is well known that the 
dimension o\ = dim T-Li is given by 

n + d-l\ fl + d-& 



d- 1 J \ d-1 

The restrictions of Y G Hi on S' d_1 are the spherical harmonics. Let oj denote the 
Lebesgue measure on and let ujd := ^(S^ 1 ) — 2n d / 2 /T(d/2). Throughout 

this section, we use { Yj (x) : 1 < j < o~i } to denote an orthonormal basis for Hi on 
S d ~ x . Thus, we have 

-/ 

Ud Js d 



Yh(x')Yj*(x')dw(x') = S h , l2 6 juj2 . (5.1) 
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Recall that in polar coordinates we have 

d- 1 



A, 



3 



1 . 
— Agd-l, 



P P- 

where Agd-i is the Laplace-Beltrami operator on the sphere S l . Since 

Y l Ax) = p l YUx'l 



the polynomials (x) satisfy the equation 



A st -iYj(x) = -l(l + d-2)Y l (x). 



(5.2) 



(5.3) 



(5.4) 



5.2. Construction of the algebra of partial differential operators. To moti- 
vate the construction, notice that if we set /3 = 4 — 1 and if we replace z 6 [0, 1] by 
p G [0, 1], where z — p 2 then for the Jacobi measure on [0, 1] considered in Section 
[2] we obtain 

(1 - z) a z p dz = 2(1 - p 2 ) a p d ~ l dp. (5.5) 

Recall that the Jacobi polynomials on B d are defined as orthogonal polynomials on 
B d with respect to measure dp a (x) = (1 — ||cc|| 2 ) a (ix, see for instance O page 38]. If 
we use polar coordinates x — px' then, up to a scaling factor, dp a (x) is a product 
of the measure in (|5.5|) on [0,1] and the surface measure dui(x') on Note 
also that the operators T>\ and V 2 defined in the beginning of Section [4] change as 
follows: 



V x = zd z = ^pd p 
V 2 = zdl + 08 + 1)0 
Moreover, in polar coordinates we have 

pd P 



(5.6) 



d 



a; • V x , 



where Va; is the gradient, while the operator d 2 



-d„ is the radial part of the 



Laplace operator A x . It is easy to see that the operators \x ■ V x and jA x satisfy 
the commutativity relation 



1 1 
-A x ,-x-V a 



= t4 
4 



which combined with (|4.3p shows that there is a natural isomorphism between the 
algebra Dp given in (I4.2[) and the associative algebra generated by \x- V x and kA x 
defined by 



Z>i = z0, 
2?2 = ^0? 



> 2 s ' V * 
08 + 1)0, 



(5.7a) 
(5.7b) 



Thus if we set 
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we see that the commutative algebra T> a ' d / 2 1;a defined in Theorem 14.21 is isomor- 
phic to a commutative algebra of partial differential operators: 



IC d (a; a) = jfl, • V, , ±A X \ : f e A d 



{a; a) 



(5.8) 



The main point now is that we can write a basis for V which diagonalizes the 
operators in IC d (a;a), using the polynomials q defined in Proposition 14.71 and the 
spherical harmonics. 

Theorem 5.1. For n G No, i G No, such that i < f and j G {1, 2, . . . , cr n _ 2 j} 
define 

QnM = Cn-^^^lkll 2 )^^ 21 ^)- (5-9) 

Then the polynomials {Qn,i,j(x)} form a basis for V d and for every f G A d (a;a) 
we have 



1 



f 



(5.10) 



Proof. The linear independence of Q n ,i,j follows easily from (I5.1[) and the fact that 
the polynomials {qf' s }i S N are linearly independent. Since 



L"/2J 



n + d-l 
d-1 



= the number of monomials of total degree n, 



we see that the polynomials Q n ,i,j form a basis for V d . It remains to show that 
they satisfy (|5.10p . Using polar coordinates we find 

Q»,M = <„^ 1; VV"- 2 ^- 2 V). 

From equations (|5.2p and f|5.4p we see that 



dl + - — -d p - \(n - 2i)(n - 2i + d - 2) 
P P 



C;^ 1; V)p"- 2i ^- 2 V)- 



The proof of (|5. 10[) now follows from Proposition 14.71 upon changing z — p 2 and 
using equations (|5.6j) . □ 



6. An explicit example 

In this section we illustrate all steps in the paper with the simplest nontrivial 
case a = k = 1. 

6.1. Krall polynomials. Let us first consider the one-dimensional case which 
leads to Krall polynomials 14 j. When a — k — 1 the polynomials q n (z) con- 
structed in Section [3] are given by the following formula 



,i,/3;»o 



(0);l,/3;a li/3 



(6.1) 



where pj l '^(z) are the Jacobi polynomials in (|2.1[) and 



(0);l,/3;a 



ao 



(n + l)(n + j9 + l) 
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For n ^ m they satisfy the orthogonality relation 



C q^{z)q^{z)z^dz + 1 q y> a °(l)q l J-' a °(l) = 0. (6.3) 
o a o{P + tj 



Since k = 1 formula g3J shows that r„ = ^ 0);1 ^ ;a ° = a + ( " +1) j"{ ;3+1) . From 

this and Remark |4. II it follows that the commutative algebra A 1 ' 13 ' 0,0 defined in (|4.7[) 
is generated by two polynomials of degrees 2 and 3. A short computation shows 
that 

A i,P;a = R[/ 2i / 3 ] ) ( 6 .4) 

where 

h{t)=t 2 + + 4ao + 4/3 + 4a /3)t (6.5a) 

/ 3 (t) = t 3 + hi + 6a„ + 6(3 + 6a /3)t 2 - —(21 + 12a + 28/3 + 12a /3 + 4£ 2 )i. 
4 16 

(6.5b) 

The algebra £) 1 '' 3;ao defined in Theorem 14. 2 1 is generated by the operators B 2 := B/ 2 
and i?3 := B/ 3 of orders 4 and 6 respectively. The operator B 2 goes back to the 
work of Krall [H] . Using the operators T>x and T>% given in (|4.ip we can write B 2 
as follows: 

B 2 (V 1 ,V 2 ) =V\- 2V 2 V\ + V\ + 2(1 + /3)X>? - 2#E> 2 £> 1 

+ (1 + 2a + 3/3 + 2a /3 + /3 2 )D 2 - 2(1 + a Q + a f3)V 2 (6 _ g) 

+ (1 + + 2a (l + /3))2?! - 4(3 + 2/3)(3 + 6/3 + 8a (l + /3)). 

16 

One can write a similar formula for B%. If we denote 

s _ 4a + 4a /3 + 2/3s + s 2 
° ~ 4(l + /3 + S ) ' 

then from formula (16. 2[) it follows easily that 



, (0);l,/3;a = 1 + P + g , (0) ; l,/3+s;ng 
Vri+s/2 1 + /3 



i.e. the functions V^+j/f' 00 ana ^ -0^°- l ' 1 '^ +s ' ao differ by a factor independent of n. 
Therefore for the polynomials q defined in (I4.26|) we find 

This combined with (I6.3[) shows that for n ^ to the polynomials ;a ° (z) satisfy 
the orthogonality relation 



2/3s 



4a (/3 



s tt) f tii' a °(z)q^r(z)ze+ s dz+ fB.^ ^wtigrw = o. 

(6.7) 
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6.2. Krall polynomials in higher dimension. Let us consider now x £ M. d and 

set (3 = I — 1. Then the algebra K, d (a,a ) defined in (I5.8[) is generated by the 
operators B2 {\x ■ V ' x , -|A X ) and B3 (J^x ■ V ' x , ^A^) which act diagonally on the 
basis of polynomials Q n .i,j described in Theorem 15.11 Let us denote by Uq the 
constant 

1 2 
U ° a (/3 + l) a d' 
Then using equations (|5.1|) . (|5.4j) . (|5.5[) and (|6.7|) we see that the polynomials Q n ,i,j 
are mutually orthogonal with respect to the inner product on V d defined by 

(/,<?>=/ f(x)g(x)dx + ^ [ f(x')g(x')Mx') 

uq • /S " 1 (6.8) 

The interesting new phenomenon in the multivariate case is the fact that, even in 
the simplest example (a = k = 1), the polynomials are orthogonal with respect to 
an inner product involving the spherical Laplacian. In this respect, the multivariate 
analogs of Krall polynomials discussed here are related to the so called Sobolev or- 
thogonal polynomials, see for instance |18j and the references therein. However the 
appearance of A S d-i in the inner product defined in (|6.8|) . which comes naturally 
from our approach, seems to be new. 

It would be interesting to find explicit orthogonality relations for the general 
multivariate polynomials Q n .i,j defined in Theorem 15. II The key step would be to 
discover an orthogonality relation similar to (|6.7j) for the polynomials <?"'f ;a given 
in Proposition ^. 71 
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